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Using the Mellin transform approach, it is shown that, in contrast with integer-order 
derivatives, the fractional-order derivative of a periodic function cannot be a function 
with the same period. The three most widely used definitions of fractional-order deriva- 
tives are taken into account, namely, the Caputo, Riemann-Liouville and Grunwald- 
(-H I Letnikov definitions. As a consequence, the non-existence of exact periodic solutions in 

a wide class of fractional-order dynamical systems is obtained. As an application, it is 
emphasized that the limit cycle observed in numerical simulations of a simple fractional- 
order neural network cannot be an exact periodic solution of the system. 
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00 ' 1. Introduction 

o . 

Although fractional calculus has a more than 300 year long history, its practical ap- 
plications are just a recent focus of interest. In the past decade, scientists and engineers 
became aware of the fact that the description of some phenomena is more accurate when 
KJi \ the fractional derivative is used. It has been recently found that many systems in intcrdis- 

j_j ■ ciplinary fields (such as acoustics, mechanics, electromagnetism, heat transfer, electrical 

' circuits, signal processing, system identification, control and robotics, chemistry, biology, 

physics, economy and finance) can be successfully described by fractional differential 
equations. 

Applications of fractional calculus and fractional-order differential equations include, 
but are not limited to modeling real world phenomena such as: dielectric relaxation 
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phenomena in polymeric materials [s^] , transport of passive tracers carried by fluid flow 
in a porous medium in groundwater hydrology |33| , viscoelastic behavior [l9| , transport 



dynamics in systems governed by anomalous diffusion [18|, |27[ , self-similar processes such 
as protein dynamics long-time memory in financial time series using fractional 
Langevin equations [Jetc. In recent years, even fractional-order models of happiness ^33 | 
and love [l7| have been developed, and they are claimed to give a better representation 
than the integer-order dynamical systems approach. 

Highly remarkable scientific books which provide the main theoretical tools for the 
qualitative analysis of fractional-order dynamical systems, and at the same time, show 
the interconnection as well as the contrast between classical differential equations and 
fractional differential equations, are |23i, 0,[3l[. 

The existence of periodic solutions is a often a desired property in dynamical sys- 
tems, constituting one of the most important research directions in the theory of dynam- 
ical systems, with applications ranging from celestial mechanics to biology and finance. 
The existence of weighted pseudo-almost periodic solutions of fractional-order differential 
equations has been investigated in 0, [l^ill^- Moreover, several results concerning the 
existence of solutions of periodic boundary problems for fractional differential equations 
have been recently reported in 0, 0, [H, 113 ■ 

In this paper, we will first show that, unlike in the case of the integer-order derivative, 
the fractional-order derivative of a periodic function cannot be a function with the same 
period. As a consequence, we obtain the non-existence of exact periodic solutions for a 
wide class of fractional-order differential systems. Finally, we exemplify the theoretical 
results in the framework of fractional-order neural network models. 



2. Preliminaries 

2.1. Fractional order derivatives 

In general, three different definitions of fractional derivatives are widely used: the 
Grunwald-Letnikov derivative, the Riemann-Liouville derivative and the Caputo deriva- 
tive. These three definitions are in general non-equivalent. However, the main advantage 
of the Caputo derivative is that it only requires initial conditions given in terms of integer- 
order derivatives, representing well-understood features of physical situations and thus 
making it more applicable to real world problems. 

Let C" (J, K) denote the space of n-times continuously differentiable functions on the 
real interval /. 

Definition 1. Let a G (0, cxd) \ N and the function g £ C"([a, 6], M). We define: 

• the Riemann-Liouville fractional- order derivative of order a of g, given by 

where n — [a] -f 1 . 

• the Caputo fractional- order derivative of order a of g, given by 

^DZ^g{t) = ^ [\t - sr-"-^g^"\s)ds, 
T(n - a) 
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where n = [a] + 1. 

• the Grunwald-Letnikov fractional- order derivative of order a of g, given by 

l{t-a)/h] 
r=0 



D2+g{t) = lim /i-" (-^y C'r js^*- 



For simplicity, when a = 0, we will drop the subscript "a + ". 

The following proposition (see 2^ 31 1) expresses the relationship between the three 
types of fractional derivatives. 

Proposition 1. Let a G (0, oo) \ N and g : (0, oo) — s- a function of class C". Then 
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D-g{t) D-g{t) D<^g{t) + ^ ^ > t^-<^ (1) 

where n = [a] + 1 . 

2.2. The Mellin transform and its properties 

In the following, we will recall the definition of the Mellin transform and a few prop- 
erties that will be helpful in the proof of our main result. We refer to 3l|, [ij, |l5| and 
the references therein, for an overview of the Mellin transform and its applications. 

Definition 2. The Mellin transform of a locally Lehesgue integrable function g : [0, oo) — >■ 

C is defined by 

poo 

M{g){z) = / g{t)e-^dt. 



The largest open vertical strip of the complex plane, of the form 

S'g = {z e C : a < < &}, 

in which the integral converges, is called the strip of analyticity ( or fundamental .strip ) 
of the Mellin transform. 

Proposition 2 (Mellin transform of tiie convolution). Considering the Mellin con- 
volution g * h of two functions g, h : [0, oo) — C defined by 



ig*h){t) = / g{ts)his)ds, 
Jo 

the following equality holds: 

M{g * h){z) ^ M{g){z) ■ M{h){l - z) 

for any z from the fundamental strip of A4{g), .such that 1 — z belongs to the fundamental 
strip of M{h). 
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Proposition 3 (Inversion of the Mellin transform). Let g : [0,oo) C be an in- 
tegrable function, such that its Mellin transform A4{g) has the fundamental strip Sg — 
{2; e C : a < 3fJ(z) < 6}. If c £ {a,b) is such that A4{g){c + is) is integrable, then the 
following equality holds: 



-1 />C+200 

— / M{g){z)t-^dz^g{t), 



almost everywhere on (0, 00). 



3. Main result 

For completeness, we first prove the following standard mathematical result: 

Lemma 1. Let n G N and T > 0. If x : (0,oo) — )• M is a non-constant T-periodic 
function of class C" on (0, oo), then for any k £ N, k < ii, the k-th order derivative x^*^^ 
is also a non-constant T-periodic function. 

Proof. Let /c G N, fc < n. If x : (0, 00) R is a non-constant T-periodic function, it 
follows that 

x{t + T)^x{t) Vie(0,cx)). (2) 
Differentiating k times in this equality, with respect to t, it follows that 

x(''\t + T)=x'-''\t) Vie (0,00), (3) 

which means that x^*^-' is also T-periodic. 

For the second part of the proof, we will proceed by reductio ad absurdum. Assuming 
that there exists fc e N, fc < n, such that x'^'^^ (t) is a constant function, we obtain that 
x{t) is a fc-th degree polynomial function 

Xit) = Cfet'' + Ck-it''-^ + ... + at + Co, 

where Ci G M, for i = 0,k. Denoting 

y{t) = x{t) - Co = Ckt'' + Ck-it''-^ + ... + cit, 

it is clear from the fact that x{t) is T-periodic, that the function y{t) is also T-periodic. 
Moreover, we note that y(0) = 0, and from T-periodicity, it follows that 

y{mT) = y{0) = Vm £ N. 

Since T > 0, this means that the fc-th degree polynomial function y{t) has an infinity of 
roots: {mT, m G N}, which is absurd. Hence, our assumption is false, and therefore, for 
any fc G N, fc < n, the function a;^'^-' (t) is non-constant. □ 

Remark 1. In the proof of the previous Lemma, when we differentiate equality ([2]) 
with respect to t, we implicitly rely on the well-known chain rule for the integer-order 
derivative: 

[/(ffW)]' = 3'W-/'(.9(i)). 
4 



This allows us to differentiate the left hand side of ([2]), leading to: 

^x{t + T)^ ^{t + T)- x'{t + T) = 1 ■ x\t + T) = x'{t + T), 
at at 

and hence: 

x'{t + T) = x{t) Vte(0,oo). 

By successive differentiations, we obtain equality 

However, it is well known that there is no simple chain rule for fractional-order deriva- 
tives (see, for example, section 2.7.3 in [sij). Hence, applying the fractional derivative 
operator with a G (0, oo) \ N, to both sides of the equation ([5]), will clearly not lead 
us to 

{D°'x){t + T) = {D°'x){t) Vie(0,oo). 

In the following, we present the main result of this paper. 

Theorem 2. Let a G (0, oo) \ N and n = [a] -I- 1. If x : {0, oo) M. is a non-constant 
T -periodic function of class C", then its Caputo derivative ^ D^x cannot be a T -periodic 
function. 

Proof. Assume that there exists a non-constant periodic function x{t) of period T > 
0, n-times continuously diffcrentiable on (0,cxd), such that its Caputo derivative is T- 
periodic as well, i.e. 

^D°'x[t)^'^ D°'x{t + T) Vt>0. 
Applying the definition of the Caputo derivative, we get: 

t rt+T 

{t - s)"~"-ia;(") {s)ds = / {t + T- s)"-"~ia;(") (s)ds, 



where n = [a] -f 1 . Making the change of variables s -\- T — s' m the integral from the 
left hand-side of the equality, and taking into account that a;(") is T-periodic as well (see 
Lemma [1]), we obtain: 

t+T pt+T 

{t + T- s)"-"-ia;(") {s)ds = / {t + T- s)"-"-^x^'''> {s)ds, 

Jo 

which is equivalent to 

{t + T- s)"-"-^a;(")(s)ds = , Vt > 0. 

Next, we make the change of variables s — T — ts' and obtain 

{t + ts)"-"-ia;(")(T - ts)tds = , V i > 

and, simplifying i"^", this is equivalent to 

(1 + s)"-"-ia;(") (T - ts)ds = , V i > 0. 
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x(")(r-it) , if u e [o,r] 

, if u > T 



Defining the function 

hiu) = 

the previous equation can be written as 

(1 + s)"-"-i/i(ts)ds = Vt>0, 



which translates into the fact that the Mellin convolution h* g of the functions h{t) and 
g{t) = (1 + t)^~°'~^ is equal to 0. Applying the Mellin transform to this equality, it 
follows that 

M{h*g){z) ^0, 

for any z where this Mellin transform is defined. 

Using the formula for the Mellin transform of the convolution given by Proposition 
[21 we get 

H{z)G{l-z) =0, 

where H = M{h) and G = M{g). This equality is true for any z G C such that z G Sh 
and 1 — z € Sg, where Sh and Sg are the fundamental strips of analiticity of the Mellin 
transforms of h and g. 

It can be easily seen that the Mellin transform of g is defined only on the strip 
Sg = {z eC -.0 < ^{z) < a - n + 1} and 

Giz)= Tgitr-^dt^ fd+t)— + 

Jo Jo r(a-n+l) 

On the other hand, evaluating the Mellin transform of h, we get 



H{z) 



nOO nl 

Jo Jo 



As is continuous and periodic, it follows that it is bounded, and we denote ||oo = 

sup Hence, 

te[o,T] 

Jo Jo 

The integral from the right hand side of the inequality is convergent if and only if 3fi(z) > 
0, and so the whole half-plane 3f?(z) > is included in Sh- 

Since H{z)G{l — z) = 0, for any z G C such that z G Sh and 1 — z e 6*9, it follows 
that 

H{z)G{l - z) = Vz e C such that n-a< 5R(z) < 1. 
Since G{1 — z) ^ for any z in the strip n — a < 3fJ(z) < 1, it follows that 

H(z) = Vz e C such that n- a< ^{z) < 1. 



Applying the inverse Mellin transform, according to Proposition [3] we obtain that the 
function h is equal to almost everywhere on (0,oo), and therefore, a;*-"-* is identically 
null (since it is continuous). Since all the derivatives of the function x, up to the order 
n, are T-periodic functions, it can be easily seen that the function x is constant, which 
contradicts the initial assumptions. The proof of the theorem is now complete. □ 

Remark 2. The proof of Theorem [5] relies on the definition of the Caputo derivative of 
order a £ (0, oo) \ N and the properties of the Mellin transform. When we talk about the 
usual integer-order derivatives, it is obvious that the derivative of a differentiable periodic 
function will be a periodic function of the same period. Therefore, Theorem [5] highlights 
a very important difference between fractional-order and integer-order derivatives. 

We emphasize that the arguments from the proof of this theorem do not hold in the 
case when a G N, i.e., when we deal with the usual integer-order derivatives. Indeed, if 
we assume that a G N, then n — [a] + 1 — a + 1 and the function g defined in the proof 
of Theorem [5] is g{t) — 1, for any t G (0, oo). Since the integral t^~^dt is divergent, 
for any z S C, it follows that the Mellin transform of the function g does not exist. 

Remark 3. In the papers [H, [s^, the authors attempt to prove the same result as the 
one presented in Theorem [21 but their arguments are flawed. One may clearly expect 
that at a certain point, the arguments provided in the proofs of the statements will fail 
for positive integer values of the order a (due to the facts described in the previous 
remarks). Nevertheless, a careful examination reveals that this is not the case in these 



papers. In fact, there is a mistake in the proof of Lemma 3 from [36|, which constitutes 
one of the pre-requirements of the main results formulated in the Theorem from fsdl and 
Theorem 1 from [35|. 

More precisely, the error comes from the fact that in the proof of the above mentioned 
Lemma 3, the areas and S~ are wrongly considered to be independent of p. However, 
since and S~ denote the area between the curve {pT — tY^~^gi{t) and the t-axis, 
located above and bellow the axis, respectively (the interval [0,T] is considered), they 
clearly depend on the values p € N. In fact, eq. (31) from fSG^ should be written as 



{pT~tf^-^Ut)dt = St{p)-S~{p). 
Therefore, based on (10) and (32) from [s^l, some simplifications lead us to: 



P 



ai-Pi Q+ 



stip) -{p- ir-^^s-{p) <o<{p- ir~^^s+{p) ~ p"^-^^s-{p), 



where at < This can be re- written as 

-ly— '^^+(p) / p y--"- 

p J ~ s~ip) -\p-ij 

Passing to the limit when p — >■ oo, we can simply obtain that 

hm m . 1 

P-^oo S- (p) 

but it is impossible to conclude that S^{p) — S~{p), for any p e N, which would lead us 
to the conclusion of Lemma 3. 
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In conclusion, it has to be underlined that the statements of Theorem 1 from [35| 
and the main Theorem from 36[ are correct, but the proofs of these statements rely on 
Lemma 3 from 36], which is incorrect. 

We emphasize that the Mellin transform approach presented in this paper for the 
proof of Theorem [2] is essentially different from the one attempted in 3^ 3^ . 



Example 1. Considering a G (0, 1) and the function x{t) = sin(t), based on Theorem 6 
from [21] we obtain 

where Ea^^ denotes the two-parameter Mittag-Leffter function defined by 



k=0 



r{ak + p) ' 



From the properties of the Mittag-Leffter function (see eg. (1.93) in [31]) it follows that 

^D" sin(i)=tl-"S2,2-a(-<'). 

For a G (0, 1), this function is not periodic (it can be verified by numerical simulations, 
see Fig. 1 ), even though it tends asymptotically to the periodic function sin + ^) • For 
a = l, we have t^-°'E2.2-a{-t'^) = E2.i{-t^) = cos(t). 




Figure 1: The functions cos {t) (blue, dashed) and t^'" E2,2-a{-t^) (red) for a = 0.5. 



The following corollary (see also Theorem 2 from [33) holds: 

Corollary 1. Let a € (0, oo) \ N and n = [a] -\- 1. If x : (0, oo) — )■ R is a non-constant 
T -periodic function of class C" on (0,oo), its Riemann-Liouville derivative ^^D'^x and 
its Grunwald-Letnikov derivative '-^^D°'x cannot be T -periodic functions. 

Proof. Based on Proposition [1] we have 



n-l 



RL 



D"x{t) =^ D"xit) + J2 ' ^ ^fc-" 



fe=0 



r(fc - a + 1) 



where n = [a] + 1. Denoting ||2:("^||oo = sup \x^^^\t)\, we have: 

te[o,T] 

i (n - a) Jo 1 (n - a + 1) 

and it is easy to see that D°'x{0^) = 0. 

If there exists k e {0, 1, n - 1} such that x^''^ (0+) ^ 0, it fohows that 

t->o+ ^' t^o+^r(fc-a+l) 

Hence, ^^D°'x{t) is unbounded in a neighborhood of 0, and it cannot be periodic. 

On the other hand, if a;('^'(0+) = 0, for any k G {0,1, — 1}, it foUows that 
^^L»"a;(t) =c D'^xit), and hence, based on Theorem^ it fohows that ^^D"x{t) cannot 
be T-periodic. 

The conclusion about the Grunwald-Letnikov derivative '-^^D°'x{t) can be easily de- 
rived from the equality ^^D°'x{t) D"x{t). □ 



4. Non-existence of periodic solutions in a class of fractional-order dynamical 
systems 

Let a = (ai, a2, ap), with ai G (0, oo) \ N for any i G {1, 2, ...,p} and m G N. We 
consider the following general class of fractional-order differential systems given in vector 
form: 

i?"x(t)-f(t,x(t),x'(t),...,x('")(t)), t>0, (4) 
where f : [0, oo) x Rp('"+i) — j. Rp is a given function and 

I?"x(t) = {D"^xi{t),D'''x2it),...,D"^Xpit)f 

where 0"% i G {l,2,...,p}, is understood as one of the fractional-order derivatives 
(Riemann-Liouville, Caputo or Grunwald-Letnikov) given in Definition 1. 

In the following, we denote n = ([ai], [a2], [c^p]) G N'' and we say that the function 
x(t) = (xi (t) , X2 (0 , • ■ • , Xp {t)Y is of class C" (/, W) if and only if Xi (t) is of class C"- (/, M) , 
for any i G {l,2,...,p}, where / denotes a real interval. 

The main results presented in the previous section lead us to the following conse- 
quence: 

Corollary 2 (The non-autonomous case). Let fc G N* and assume that the function 
f is T-periodic with respect to its first argument. Then there are no non-constant kT- 
periodic solutions of class C" of system 

Proof. Assuming that there is a non-constant fcT-periodic solution x(t) of class C" of 
system from the T-periodicity of the function f with respect to its first argument, it 
follows that _D"x(t) is fcT-periodic as well, which contradicts Theorem [2] or Corollary [T] 
□ 
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In a very similar manner, the following result can be proved in the autonomous case: 

Corollary 3 (The autonomous case). // the function f is constant with respect to 
its first argument (i.e., the system is autonomous), then there are no non-constant 
periodic solutions of class C" of system 

Remark 4. Corollaries [2] and [3] hold even if some (but not all) of the differentiation 
orders ai,a2,...,ap from system ^ are positive integers, i.e., some (but not all) of the 
fractional derivatives are replaced by integer-order derivatives. 

Remark 5. Even though, based on Corollary [31 exact periodic solutions do not exist 
in autonomous fractional-order systems, oscillatory behavior (limit cycles) has been ob- 
served by numerical simulations in many systems such as: a fractional-order Van der Pol 
system fractional-order Chua and Chen's systems , [ll[ , a fractional-order Rossler 
system '38] and a fractional-order financial system 

Remark 6. In the past few years, important existence results for solutions of periodic 
boundary value problems involving fractional differential equations have been obtained 0, 
I37I I , but the theoretical setting is essentially different from the one presented in Corollaries 
[Ilandll 

Remark 7. We emphasize that it is possible to obtain exact periodic solutions in im- 
pulsive fractional-order dynamical systems, by choosing the correct impulses at the right 
moments of time. For example, we may consider the impulsive dynamical system 

^7^"x(t) =/(t,x(i)), t>0,t^tk, 

Ax|^^^^ =4(x), fcGZ; (5) 

x(0) = xo, 

where a £ (0,1),/: [0, 00) x M" — > R" is a given function, xq G M", the sequence of times 
{tk)kez+ is strictly increasing, to = 0, Ax\t=t^ = x(i+) - x(t^), where x(i+) and x(i^) 
represent the right and left limits of x(t) at t — tk, and Ik are the impulsive operators. 
We consider that the following assumptions are satisfied: 

• The function f{t, •) is T-periodic. 

• Let p e Z+ such that [0,T]n(tfe)feGZ+ — {to,ti,t2, ■■■,tp}- We assume: t^+p — tk+T, 
for any k G Z+. 

Defining the impulsive operators by 

4(x) = --i- / ' {tk - s)"- x(s))ds, Vfc e z; 
r(a) 

ensures that the solution of ([5]) is T-periodic. Theoretical details and several applications 
will be given in a future paper. 
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5. An application to fractional-order neural networks 



The main advantage of fractional-order mathematical models in comparison with clas- 
sical integer-order models is that fractional derivatives provide an excellent tool for the 
description of memory and hereditary properties of various processes. In fact, fractional- 
order systems have infinite memory. Taking into account these facts, it is easy to see that 
the incorporation of a memory term (in the form of a fractional derivative or integral) 
into a neural network model is an extremely important improvement. 

Based on this idea, the common capacitor from the continuous-time integer-order 
Hopfield neural network is replaced by a generalized capacitor, called fractance, giving 
birth to the so-called fractional-order Hopfield neural network model. 

The fractional-order formulation of artificial neural network models is also justified 
by research results concerning biological neurons. In the recent paper [2^ it has been 
pointed out that fractional differentiation provides neurons with a fundamental and gen- 
eral computation ability that can contribute to efficient information processing, stimulus 
anticipation and frequency-independent phase shifts of oscillatory neuronal firing, em- 
phasizing once again the utility of developing and studying fractional-order mathematical 
models of neural network dynamics. 

The analysis of fractional-order artificial neural networks is a very recent and promis- 
ing research topic. The first papers discussing fractional-order neural network models 
0, 26, [2^, Q report on results of numerical simulations, especially on the numerical 
evidence of limit cycles and chaotic phenomena. 

We consider a simple fractional-order neural network of two neurons, given by the 
following autonomous system of fractional-order differential equations 



^D"xi{t) = -xi{t) + 2tanh(a;i(<)) - 0.5 tanh(a;2(i)) 
^D°'x2{t) = -X2it) + tanh(.xi(i)) + 2 tanh(.T2(t)) 



(6) 



where the fractional order is a = 0.5. 




Figure 2: The trajectory of system with the initial condition {xo,yo) = (0.1,0.1) converges to an 
asymptotically stable limit cycle. 



Numerical simulations suggest that there exists an attracting limit cycle for system 
dl]) (see Figures 2 and 3). However, based on Corollary [31 the attracting limit cycle is not 
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Figure 3: 
(0.1,0.1). 



The solution {xi{t), X2{t)) of system ^ versus time, with the initial condition {xo,yo) 



an exact periodic solution of system 
system 



We refer to [22] for additional results concerning 



6. Conclusions and directions for future research 

In this paper, it has been shown that the fractional-order derivative of a periodic 
function cannot be a periodic function of the same period, underlining a very important 
difference between integer and fractional-order derivatives and explaining the absence of 
periodic solutions in a class of fractional-order dynamical systems. Since this remarkable 
property holds for three widely used types of fractional derivatives (Riemann-Liouville, 
Caputo and Grunwald-Letnikov) , it may be interesting to investigate if it is also true for 
other known definitions of fractional derivatives, or if it is possible to give a new formal 
definition of a fractional-order derivative which preserves periodicity. For example, in 
section 3.1 of a generalized definition was presented for fractional derivatives of 
periodic functions, using Fourier expansions, and several results concerning periodicity 
were also give n. Furthermore, two kinds of Weyl fractional derivatives defined for periodic 
functions [20|, namely the Weyl-Liouville and Weyl-Marchaud fractional derivatives, are 
also worth mentioning. 
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